Elucidating the Vacuum Structure of the Aoki Phase 



V. Azcoiti 1 , G. Di Carlo", E. Follana 1 and A. Vaquero" 1 . 
October 5, 2012 

Abstract 

In this paper, we discuss the vacuum structure of QCD with two flavors 
of Wilson fermions, inside the Aoki phase. We provide numerical evidence, 
coming from HMC simulations in 4 4 , 6 4 and 8 4 lattices, supporting a vac- 
uum structure for this model at strong coupling more complex than the one 
assumed in the standard wisdom, with new vacua where the expectation value 
of iyjsy can take non-zero values, and which can not be connected with the 
Aoki vacua by parity-flavour symmetry transformations. 



^ 1 Introduction 

43' 



Two and three colors QCD with unimproved Wilson fermions started to be simu- 
lated in the early 80's (HQ]. The complexity of the phase structure of this model 
was known from the very beginning, and indeed the existence of a phase at fi- 
nite lattice spacing, a, with spontaneous parity and flavor symmetry breaking, was 
conjectured for this model by Aoki in the middle 80's Jl, H]. From that time on, 
much work has been done in order to confirm this pattern of spontaneous symme- 
try breaking, and to establish a quantitative phase diagram for lattice QCD with 
Wilson fermions. References §5- 26] are an incomplete list of the work done along 



these years. 

Aoki's conjecture has been supported not only by numerical simulations but 
also by theoretical studies based on the linear sigma model I6J] and on app lying 



chiral perturbation theory (%PT) to the continuum effective lagrangian CI 111 . The 
latter analysis predicts, near the continuum limit, two possible scenarios, depend- 
ing on the sign of an unknown low-energy coefficient. In one scenario, flavor and 
parity are spontaneously broken, and there is an Aoki phase with a nonzero value 
only for the ityYsXsy condensate. In the other one (the first-order scenario) there is 
no spontaneous symmetry breaking. 
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This standard picture for the Aoki phase was questioned three years ago by 
three of us in 112 ill , where we conjectured on the appearance of new vacua in the 
Aoki phase, which can be characterized by a non-vanishing vacuum expectation 
value of the flavor singlet pseudoscalar condensate /\j/y5\|/, and which can not be 
connected with the Aoki vacua by parity-flavor symmetry transformations. How- 
ever Sharpe pointed out in 02211 that our results, based on an analysis of the Prob- 
ability Distribution Function (PDF) of fermion bilinears, could still be reconciled 
with the standard scenario, and that otherwise we would question also the validity 
of the %PT analysis. 

For the last few years we have performed an extensive research on the vacuum 
structure of the Aoki phase 123142511 . in order to find out if our alternative vacuum 
structure, derived from the PDF analysis, was realized or not. The purpose of this 
paper is to clarify these issues by reporting the results of our investigations which, 
as will be shown along this article, provide evidence of a more complex vacuum 
structure in the Aoki phase of two-flavor QCD, as conjectured in 112 IP . 

The outline of the paper is as follows. Section 2 summarizes the standard pic- 
ture on the phase diagram of QCD with two flavors of unimproved Wilson fermions 
and our alternative scenario, derives some interesting expressions to be compared 
with continuum results, and analyzes the relevance of a discrete symmetry P ', com- 
position of parity and discrete flavor transformations, which was introduced by 
Sharpe and Singleton 111 ill as an attempt to reconcile the Vafa-Witten theorem on 
the impossibility to break spontaneously parity in a vector-like theory, with the 
existence of the Aoki phase. Section 3 contains technical data of our HMC simu- 
lations of QCD with two flavors of Wilson fermions, inside and outside the Aoki 
phase, performed without external sources and also with a twisted mass term in the 
action. The results of our analysis of the PDF of the operator Q, an order param- 
eter for the P' symmetry, are reported in Section 4. This Section contains what is, 
in our opinion, the strongest indication on the existence of a vacuum structure in 
the Aoki phase more complex than the standard one. In Section 5 we show how 

we were able to perform a direct measurement of ^(/\j/y5\|/) 2 ^) in the Aoki phase, 
which produced a non-zero value for this operator in the 6 4 lattice, of the same or- 
der of magnitude as ( (/\j/„Y5\|/ u ) 2 ) and ( (/\j/y5X3\|/) 2 ) (both are non-vanishing in 



the Aoki phase in the standard picture). Finally Section 6 contains our conclusions. 



2 The Aoki phase for the two flavor model 

The fermionic Euclidean action of QCD with two degenerate Wilson flavors is 

S F = %W(K)y u +\\r d W(K)\\f d , (1) 

where W (k) is the Dirac -Wilson operator, and K is the hopping parameter, related 
to the adimensional bare fermion mass niQ by K = 1/(8 + 2mo). The standard 
wisdom on the phase diagram of this model in the gauge coupling P,K plane is the 
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one shown in figure Q] The two different regions observed in this phase diagram, 
A and B, can be characterized as follows: in region A parity and flavor symmetries 
are realized in the vacuum, which is supposed to be unique. The continuum limit 
should be obtained by taking the g 2 — > 0, K — > 1/8 limit from within region A, 
which we will call the QCD region. In region B, on the contrary, parity and flavor 
symmetries are spontaneously broken, there are many degenerate vacua connected 
by parity-flavor transformations in this region, and the Gibbs state is therefore made 
up from many degenerate vacuum states. In what follows we will call region B the 
Aoki region. 



CO 
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Figure 1: Aoki (B) and physical (A) region in the ((3,k) plane. Adapted from ifnll 
by courtesy of the authors. 



2.1 Standard Wisdom 

The standard order parameters to distinguish the Aoki region from the QCD region 
are z\jjry5T/V|/ and rij)Y5\|/, with Xj the three Pauli matrices. For the standard election, 
j = 3, they can be written as a function of the up and down quark fields as follows 

iWW = i¥«Y5V«-*'W5Yd, 
«W5¥ = «V«Y5¥« + Wsfii- 
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The standard wisdom for the Aoki phase can be summarized by the following 
two equations |0] 



(z'WaV) + 0, 

These equations should hold in the vacuum selected when we add a twisted 
mass term to the euclidean Wilson action (Q]), 



S m , = m f jYY 5 t 3 \|/, (4) 

which breaks both parity and flavor symmetries, and take the vanishing twisted 
mass limit once the infinite degrees of freedom limit has been performed. The 
first of the two condensates iyjsTiy breaks both parity and flavor symmetries. 
The non-vanishing vacuum expectation value of this condensate signals the spon- 
taneous breaking of the SU (2) flavor symmetry down to U(l), with two Goldstone 
charged pions. Notwithstanding parity is spontaneously broken if the first equa- 
tion in © holds, the vacuum expectation value of the flavor diagonal condensate 
(/\jjry5\|/) vanishes. Indeed it is also an order parameter for a discrete symmetry P ', 
composition of parity and discrete flavor rotations, which is assumed to be realized 

iH. 

The Aoki phase can then be characterized by the presence of two massless 
charged pions and a massive neutral pion which becomes massless just at the crit- 
ical line separating this phase from the physical one. In addition the r\ -meson, the 
mass of which can be extracted from the long distance behavior of the two-point 
correlation function of the flavor singlet pseudoscalar operator /\j/y5\|/, stays mas- 
sive even at the critical line as a consequence of the U(1)a axial anomaly. 



2.2 Our alternative Scenario 

Three years ago this standard scenario was revised in l2lll with the help of the 
probability distribution function (PDF) of the parity-flavor fermion bilinear order 
parameters I27l l28ll . We found that if the spectral density pu(k,K) of the overlap 
hamiltonian Y5W(k), or Hermitian Dirac -Wilson operator, in a fixed background 
gauge field U were not symmetric in X in the thermodynamic limit for the relevant 
gauge configurations, Hermiticity of WfuYstyu + i^djs^d would be violated at finite 
p\ Assuming that the Aoki phase ends at finite p\ this would imply the loss of 
any physical interpretation of this phase in terms of particle excitations. If on 
the contrary the spectral density pu(k,K.) of the Hermitian Dirac -Wilson operator 
in a fixed background gauge field U is symmetric in X, we argued in 112 ill on the 
appearance of other phases, which can be characterized by a non- vanishing vacuum 
expectation value of z\j/„Y5\|/ u + i^d^s^d, and which can not be connected with the 
Aoki vacua by parity-flavor symmetry transformations. We summarize here the 
main steps of 12 ill . 
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The key point was to study the Fourier transform of the PDF of the fermion 
bilinear order parameters, P{q), defined as 

= f[dU][da?d\f]exp{-SYM + ?W(k) ¥ + ffd*x 0{x)} 
[q) f[dU][d\\rd\\r]exv{-Sy M + \\rW(K)V} 

where O is any fermion bilinear operator 

0(x)=\\r(x)0\\r(x), (6) 

with O a matrix with Dirac, color and flavor indices, and V is the number of lattice 
points. Integrating out the fermion fields one gets 

= f[dU]e- s ^det(W(K) + ^d) 
(q) f[dU]e- s rMdetW(K) ' 1 ' 

which can also be expressed as the following mean value 

/ det(W(K) + |Q) \ 
P{q) = \ det^(K) /' (8) 

computed in the effective gauge theory with the integration measure 

[dU]e~ SYM detW(K). 

The q-derivatives of P(q) give us the moments of the PDF of our fermion bilinear 
d n P(q) 



= i 

q=0 



fijd'xo^yy (9) 



dq" 

For the case we are interested in, and if we call Pj(q),Po(q) the PDF of 



j (10) 
c = iyysV- 



In momentum space, we have 

^2 



p Y (q) =m-m = (u[-^ 



^-(n(vl- +1 ) 



3 ' 1 en) 



where Xj are the real eigenvalues of the one flavor Hermitian Dirac -Wilson operator 
W(K)=y 5 W(K). 

Since co,c 7 - are order parameters for the symmetries of the lattice action and 
we are not selecting a particular vacuum state, their first moment will vanish al- 
ways, independently of the realization of the symmetries. The first non-vanishing 
moment, if the symmetry is spontaneously broken, will be the second: 
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(12) 



In the QCD region flavor symmetry is realized. The PDF of C3 will be then 
8(c 3 ) and (cj) = 0. We get then 




(cl) = 

which should vanish since parity is also realized in this region. Furthermore a 
negative value of (eg) would violate Hermiticity of nj/y5\|/. 

In the Aoki region |0] there are vacuum states in which the condensate C3 (TTQb 
takes a non- vanishing vacuum expectation value. This implies that the PDF P(cj) 
will not be 8(c3) and therefore (c^) (fT2l) will not vanish. Indeed expression (fT2l for 
(c\) seems to be consistent with the Banks and Casher formula |29j] which relates 
the spectral density of the Hermitian Dirac -Wilson operator at the origin with the 
vacuum expectation value of c 3 

If, on the other hand, (/\|/y5\|/) = in one of the Aoki vacua, as conjectured in 
01, (j'WsV) = in all the other vacua that are connected with the standard Aoki 
vacuum by a parity-flavor transformation, since £\j/Ys\|/ is invariant under flavor 
transformations and changes sign under parity. Therefore if we assume that these 
are all the degenerate vacua, we should conclude that P(co) = 5(co) and (c" } ) = 0, 
which would imply an infinite tower of non trivial relations, one for each even 
moment of the PDF. We write here the simplest one, 



(13) 

A possible scenario which was discussed in Il21ll is the one that emerges if 
we assume a symmetric distribution for the eigenvalues of the hermitian Wilson 




operator. As discussed in 02111 this assumption is necessary in order to preserve 



hermiticity. If we take this assumption in the most naive way we get 




te> -<<«>= 4 TJL^T >=0, (14) 



and therefore we should conclude that new vacua characterized by a non vanishing 
vacuum expectation value of the singlet-flavor pseudoscalar operator /\j/y5\|/ must 
exist besides the Aoki vacua. Nevertheless as Sharpe noticed in Jiill . sub-leading 
contributions to the spectral density may affect equation (fT4l in the e-regime in 
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such a way that, not only (cq), but every even moment of z'\j/y5\|/ would vanish, 
restoring the standard Aoki picture. 

As we will show in this article equation (fl4l) turns out not to be realistic in 
the Aoki phase, as follows from numerical simulations. However the thesis of 
Sharpe, although possible and inspired in the absence of our conjectured new vacua 
in the chiral effective lagrangian approach, would enforce as discussed before an 
infinite series of sum rules, similar to those found by Leutwyler and Smilga in the 
continuum H30TI - The main purpose of this paper is to clarify these issues. 



2.3 Some interesting relations 

As previously discussed flavor and parity symmetries should be realized in the 
vacuum of the physical phase. This phase can therefore be characterized, in what 
concerns its low energy spectrum, by the existence of three degenerate massive pi- 
ons, which become massless at the critical line, and the r\ -meson, which is massive 
even at the critical line because of the axial anomaly. When we cross the critical 
line and enter into the Aoki phase, the neutral pion becomes massive whereas the 
charged pions are massless since they are the two Goldstone bosons associated to 
the spontaneous breaking of the SU (2) flavor group to a U(l) group. 

The susceptibilities of the neutral pion % n o and the eta-meson %^ in the physical 
phase are the corresponding integrated two-point correlation functions 




lit* =£(irowM *W3Y(o)) : , 

n 2 1 v 

Xti = £ («PYs¥ (*) 'W5 V (0) )□ 




(=1 



(15) 
(16) 



The rightmost hand sides of these equations are just the second moments of the 
PDF of ixj/YstW and z'\j/y5\|/ (see equations (fT2l ) multiplied by the corresponding 
volume factor. Therefore equations (fl5T ). and ([Tot give us the following relation 
between % n a and the trace of the inverse one-flavor hermitian Dirac -Wilson 
operator: 



Xr ] =Xno--(Tr 2 (y 5 W- l (K))), 



(17) 



to be compared with the well known relation between the eta, pion and topological 
susceptibilities 



XT| — %7t +4 r. , 



(18) 



which holds in the continuum and also in the Ginsparg-Wilson regularization at 
finite lattice spacing. 



'There are no disconnected contributions in this phase. 
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The first interesting observation that follows from equations (fTTT ). and (TT81 is 
the suggestive relation between the topological susceptibility %j and the trace of 
the inverse hermitian Dirac -Wilson operator, 

g = -i(r^( Y5 w- 1 (K))>, (19) 

which should hold near the continuum limit. This suggests also the following 
relation between the trace of the inverse hermitian Dirac-Wilson operator, quark 



mass m, and the density of topological charge q = Q/V Oil. 13211 



i = Irr(y 5 W- 1 (K)), (20) 
m V 

A second interesting observation that follows from (fPTl) concerns the behavior 

of 



l<rr 2 (y 5 lV- 1 (K))> = l 



v 1 



(21) 



when we approach the critical line from the physical phase. Indeed in the physical 
phase this quantity should be finite because the pion and eta susceptibilities are 
finite. However when approaching the critical line it should diverge in order to 
compensate the divergence of the pion susceptibility keeping %^ finite in (flTl) . in 
deep similarity with the divergence of the topological susceptibility divided by the 
square quark mass in the continuum and Ginsparg-Wilson regularization in the 
chiral limit. The origin of this divergence lies in the accumulation of eigenvalues 
of the hermitian Dirac-Wilson operator near de origin. 

The divergence of (I2TT) at the critical line is slower than V in such a way as 
to keep a parity-flavor symmetric vacuum, as corresponds to a second order phase 
transition. On the other hand if the standard scenario for the Aoki phase is realized, 
the higher accumulation of eigenvalues of the hermitian Dirac-Wilson operator 
at the origin in this phase should be enough to give a finite contribution to (c?) 
(second equation in (fT2l). Indeed if the Aoki vacuum selected with a twisted mass 
term plus those obtained by parity-flavor symmetry transformations are the only 



degenerate vacua in this phase, the following relation should hold j22l 12711 

( c 3> = y, (22) 

where £ is the mean value of c$ in the vacuum selected by the twisted mass term. 
In Section 5 we will report the results of a check of equation (l22l . 

The standard scenario for the Aoki phase requires also that the singlet flavor 
pseudoscalar operator /\j/y5\|/ takes a vanishing vacuum expectation value. Taking 
into account the first equation in (fT2l . this requirement implies that A (Tr 2 (y$ W~ l (t 
in (|2TI ) should diverge as V, to compensate the first contribution to (eg) in (fT2l) . 
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2.4 Spontaneous Symmetry Breaking of P'? 



Years ago Sharpe and Singleton suggested 111 ill that a symmetry P', composition of 
parity and flavor transformations, associated to a discrete subgroup of the parity- 
flavor continuous group, and acting on pure gluonic operators as parity, could still 
be realized in the Aoki phase, notwithstanding parity and flavor are spontaneously 
broken in this phase. Their motivation for such a proposal was the attempt to rec- 



oncile the Vafa-Witten theorem [33] on the impossibility to break spontaneously 
parity in a vector-like theory, with the existence of the Aoki phase. Their point was 
that Vafa-Witten theorem, although it does not apply to fermionic order parame- 
ters, could still apply to pure gluonic operators. Then if the realization of P' does 
not require a vanishing Aoki condensate, and P' acts on pure gluonic operators as 
parity, the existence of the Aoki phase besides the realization of P' in the vacuum 
would not be in conflict with the Vafa-Witten theorem for pure gluonic operators. 
Furthermore the realization of P' would be quite useful since it would give a sim- 
ple explanation for the tower of sum rules (fT3l) required in the standard scenario. 
Indeed the flavor-singlet pseudoscalar condensate /\j/y5\|/ is an order parameter for 
P'. 

Although the main motivation for the introduction of the P' symmetry, the real- 
ization of Vafa-Witten theorem for pure gluonic operators, has become less relevant 
on the light of later works on this subject 13411 - 13711 . the implications on the standard 
scenario for the Aoki phase of the realization of this symmetry are still relevant. 

A possible election for P', the realization of which would not be in conflict with 
non vanishing condensates (/\j/y5i ; \|/) for j = 1,2,3 is the composition of parity 
with the Z4 subgroup of the SU (2) flavor group generated by 

-^=(Xi+X 2 +T 3 ). 

In this context let us consider the operator 

Q = ±Tr(y 5 W- l (K)) = ±tv> (23) 

which is a non local order parameter for P'. This operator, although non-local, 
is not singular since exact zero modes have vanishing integration measure in the 
Wilson formulation. We will assume in what follows that Q is an intensive op- 
erator. Notice that any local operator should be intensive but any intensive op- 
erator is not necessarily local. Our assumption is based in the following argu- 
ment: (Q) = (L. Y ^(X)Y5y(X)), W{x)Y5W{x) being a local fermionic operator. 
Our prejudices tell us that ^7 LxVWYsVW' as anv intensive operator, should be 
self-averaging. In other words, if the system size is very large, a single thermalized 
gauge configuration should be enough to get ^ (L. V ¥( X )Y5V( X ))> an d this seems 
indeed to be the case in the physical phase, where the PDF of Q approaches a delta, 
as shown in Section 4. Furthermore equations (flTT ) and (fT8T ) strongly suggest that, 
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near the continuum limit, Q would be essentially the density of topological charge, 
which indeed is an intensive operator. 

Hence were the P' symmetry realized in the vacuum, the PDF of (l23l would 
be, according to our assumption, a 8(x) distribution. But the second moment of 
this PDF, 



v 1 

needs to be finite in the standard scenario for the Aoki phase in order to realize 
the first sum rule (PT3T ). Therefore the standard scenario requires the spontaneous 
symmetry breaking of P', but then there are no symmetry reasons to justify this 
tower of sum rules. In Section 4 we will report numerical results for the PDF of 
\, both in the Aoki and in the physical phase. 



3 The simulations 

In order to determine the phase structure and the properties of the vacuum, we 
need to measure the PDF's of the parity-flavor order parameters described in the 
previous sections, including the Q operator (l23l) . So we decided to carry out HMC 
simulations of QCD with two flavors of Wilson fermions, inside and outside the 
Aoki phase, and mainly without external sources, although some simulations with 
a twisted mass term in the action have also been performed. 
We remove the external sources because of several facts: 



1. The PDF formalism requires the removal of any external sources in the ac- 
tion. We just perform one long simulation in the Gibbs state. 

2. The addition of a twisted mass external source, as has been done in past 
simulations of the Aoki phase, automatically selects the vacuum where the 



standard properties of the Aoki phase are verified 12 111 . This point could 
explain why no one ever saw a new Aoki phase like the one we are proposing, 
since all the past simulations performed with two flavors of Wilson fermions 
inside the Aoki phase were done with an external twisted mass term. 

We could try to select one of these hypothetical new vacua by adding an ex- 
ternal source proportional to z\j/y5\|/, but this introduces a severe sign problem 
in the simulations. 



The Aoki phase without external sources is very hard to simulate. Inside the 
Aoki phase there are exactly massless pions and quasi-zero modes, rendering the 
condition number of the Wilson Dirac operator quite high. A standard solver will 
not be enough to invert the Dirac operator in a HMC simulation. Fortunately, in 
recent years a number of efficient solvers have appeared, and by using a DD-HMC 
113 811 we could invert the Dirac operator at a reasonable speed with the resources 
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available to us, i.e., the clusters of the department of theoretical physics of the 
University of Zaragoza, comprising 160 cores and 280 Gb of available memory, 
interconnected by a gigabit network. The simulations were parallelized for 4, 8 
or 12 cores using openMP, and we did several simultaneous runs. The iteration 
count of the solver ranged from a handful outside the Aoki phase, to a few hundred 
inside the Aoki phase in the largest volume 8 4 and without external sources. For 
volumes higher than 8 4 , the inversion time became prohibitive for our computing 
resources and our time-frame, and we should look for new ways of simulating the 
Aoki phase. 

Moreover, there is an additional problem with these simulations: since quasi- 
zero modes appear inside the Aoki phase, the eigenvalues sometimes attempt to 
cross the origin, and they would do so, were it not for the fact that the crossing 
of eigenvalues is forbidden by the HMC dynamics: the potential energy in the 
molecular dynamics step diverges at the origin, and there is an infinite repulsion 
that prevents the eigenvalues from crossing. 

In order to solve this problem we classified our simulations according to its sec- 
tor number, i.e., the number of 'crossed' eigenvalues of the hermitian Dirac -Wilson 
operator they had: beginning from a completely symmetric state (same number of 
positive and negative eigenvalues), the number of the eigenvalues that should cross 
to reach the desired state; and we performed simulations with a different number 
of crossed eigenvalues for each volume. Now this does not completely solve the 
problem, since we still don't know the weight of each sector within the partition 
function. The only way to simulate dynamically all the sectors is to add a twisted 
mass term, but as explained above we are primarily interested in the results with- 
out external sources. Fortunately the weights of the different sectors evolve very 
slowly with the twisted mass, and we can then extrapolate the data to vanishing 
twisted mass. 

Indeed in Table Q] we report the weights of the different sectors as a function 
of the volume and twisted mass m t . The values marked with * were obtained from 
reweighting and extrapolation of the configurations generated at m t = 0.01 in the 
case of the 4 4 , and m t = 0.05 in the case of the 6 4 (the m t = 0.01 simulation for 
the 6 4 came later, and has less statistics), Qch stands for quenched, whereas the 
numbers marked as MFA come from a MFA (Microcanonical Fermion Average) 
inspired approach J39I1 . and were obtained by diagonalizing 4xl0 6 gauge configu- 
rations. As can be seen in this table the weights show a very mild m f -dependence, 
although the 4 results at m t = obtained from MFA simulations are at two stan- 
dard deviations from the extrapolated results. It is not clear to us if this discrepancy 
has a statistical origin or reflects a discontinuity of the weights of different sectors 
at m t = 0. However if the latter were the actual case, the standard picture for the 
Aoki phase could not be realized since it requires continuity in the sample of rel- 
evant gauge configurations at m t = 0. Hence we will assume continuity of the 

2 Both methods agreed for our data, except for the higher sectors (2,3) of the 6 4 , where the low 
statistics prevented reweighting from working properly 
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weights, and in the following sections we will use Table [Q to reconstruct the PDF 
for the interesting bilinears by summing the weighted PDF's of all sectors for a 
given volume. The details of the runs can be checked in Table [2] 



Vol 


m t 


Sector 


Sector 1 


Sector 2 


Sector 3+ 


4 4 


MFA 


80% ± 6% 


20% ±6% 


0% 


0% 


4 4 


0.0* 


67.6% ±1.3% 


32.2% ±1.4% 


0.11% ±0.07% 


0% 


4 4 


0.01 


68.1% ±1.0% 


3 1.6% ±1.0% 


0.24% ±0.08% 


0% 


4 4 


0.05 


70.4% ±0.3% 


29.4% ±0.3% 


0.18% ±0.07% 


0% 


4 4 


0.10 


69.3% ±0.9% 


30.4% ±0.9% 


0.28% ±0.04% 


0% 


4 4 


0.30 


70.3% ±0.7% 


29.6% ±0.6% 


0.14% ±0.05% 


0% 


4 4 


0.50 


70.0% ±0.5% 


29.8% ±0.5% 


0.20% ±0.06% 


0% 


4 4 


Qch 


7 1.4% ±0.3% 


28.5% ±0.2% 


0.07% ±0.02% 


0% 


6 4 


0.0* 


36% ±4% 


59.2% ±2.3% 


4.77% ±0.45% 


0.12% ±0.05% 


6 4 


0.01 


5 1.8% ±0.6% 


44.7% ±0.7% 


3.3% ±0.3% 


0% 


6 4 


0.05 


49.9% ±0.8% 


45.7% ±0.4% 


4.3% ±0.5% 


0.08% ±0.02% 


6 4 


0.10 


5 1.4% ±0.6% 


44.7% ±0.6% 


3.86% ±0.07% 


0.06% ±0.02% 


6 4 


0.30 


52.7% ±1.2% 


44.3% ±1.0% 


2.96% ±0.20% 


0.06% ±0.06% 


6 4 


0.50 


53.1% ±0.6% 


44.0% ±0.5% 


2.82% ±0.18% 


0.06% ±0.02% 


6 4 


Qch 


55.0% ±1.6% 


43.0% ±1.9% 


2.0% ±0.4% 


0% 


8 4 


0.05 


36.0% ±0.6% 


48.2% ±0.7% 


13.9% ±0.4% 


2.06% ±0.20% 


8 4 


0.10 


36.1% ±0.6% 


47.0% ±0.4% 


14.6% ±0.3% 


2.22% ±0.28% 


8 4 


0.30 


38.2% ±0.4% 


46.8% ±0.4% 


13.33%±0.12% 


1.6% ±0.4% 


8 4 


0.50 


39.0% ±0.3% 


48.5% ±0.3% 


11. 33% ±0.25% 


1.18%±0.11% 



Table 1: Weights of the different sectors as a function of the volume and m t . The 
values marked with * were obtained from reweighting and extrapolation of the 
configurations generated at m t = 0.01 for the 4 4 and m, = 0.05 for the 6 4 . 

The high variability of the acceptance ratio is due to the hard task of fine-tuning 
the solver parameters and the trajectory length to work properly for each case. As 
the table shows, certain simulations required a very small time-step for the HMC 
to work; in those simulations the eigenvalues were trying to cross the origin, gen- 
erating high forces in the molecular dynamics, and thence creating a high rejection 
rate, so this small time-step was strictly necessary to obtain reasonable data. To 
enhance acceptance, we introduced the replay trick in the hardest calculations, i.e., 
those inside the Aoki phase and without a twisted mass term. During replays, the 
stepsize was halved whereas the trajectory length was kept constant. 

In order to calculate the PDF's we diagonalized the hermitian Dirac -Wilson 
operator on all configurations generated and found the eigenvalues. The diago- 
nalizations were carried out by using a parallelized Lanczos algorithm, combined 
with a Sturm bisection. The algorithm was checked heavily against the LAPACK 
library before starting production to ensure that our results were correct, but our 
algorithm was much faster than the LAPACK standard diagonalization procedure 
for hermitian matrices, and used a small fraction of the memory, because in our al- 
gorithm the fermionic matrix was generated on-the-fly. The diagonalization times 
ranged from around 100 seconds for the 4 4 to a few hours (close to four) in the case 
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Vol 


R 
r 


K 




Sector 


Confs 


Acc (%) 


HMC Step 




4 4 


2.00 


0.25 


0.00 


o 


5000 


99% (72%) 


3 125e-03 


0.1 


4 4 


2.00 


0.25 


0.00 


1 


2473 


99% (73%) 


2.5e-03 


0.1 


6 4 


2.00 


0.25 


0.00 


o 


676 


92% (60%) 


6 25e-04 


0.04 


6 4 


2.00 


0.25 


0.00 


1 


2000 


88% (56%) 


6 25e-04 


0.04 


4 4 


4.00 


0.18 


0.00 


All* 


10000 


92% 


0.1 


1.0 


6 4 


4.00 


0.18 


0.00 


All* 


5000 


84% 


0.1 


1.0 


8 4 


4.00 


0.18 


0.00 


All* 


5000 


74% 


0.1 


1.0 


4 4 


2.00 


0.25 


0.01 


All 


5001 


83% 


0.0133 


0.2 


4 4 


2.00 


0.25 


0.05 


All 


5001 


86% 


0.05 


0.4 


4 4 


2.00 


0.25 


0.10 


All 


5664 


92% 


0.025 


0.5 


4 4 


2.00 


0.25 


0.30 


All 


4999 


96% 


0.1 


0.5 


4 4 


2.00 


0.25 


0.50 


All 


10000 


99% 


0.025 


0.5 


6 4 


2.00 


0.25 


0.01 


All 


2520 


99% (89%) 


0.01 


0.4 


6 4 


2.00 


0.25 


0.05 


All 


5000 


70% 


0.05 


0.4 


6 4 


2.00 


0.25 


0.10 


All 


5000 


95% 


0.05 


0.4 


6 4 


2.00 


0.25 


0.30 


All 


5001 


90% 


0.1 


0.5 


6 4 


2.00 


0.25 


0.50 


All 


5001 


99% 


0.05 


0.4 


8 4 


2.00 


0.25 


0.05 


All 


5001 


77% 


0.04 


0.2 


8 4 


2.00 


0.25 


0.10 


All 


5001 


84% 


0.05 


0.5 


8 4 


2.00 


0.25 


0.30 


All 


5001 


96% 


0.05 


0.4 


8 4 


2.00 


0.25 


0.50 


All 


10000 


98% 


0.05 


0.5 



Table 2: Details of the simulations that generated our configurations. The word 
Sector refers to the eigenvalue sector where the simulation was performed, as ex- 
plained above. The * means that only sector contributes in this case, so all sectors 
were taken into account. Also, Confs refers to the number of configurations saved 
(we saved one configuration for each two generated, so as to reduce the autocor- 
relations), Acc is the acceptance ratio of the simulation (the value in parenthesis 
marks the acceptance without the replay trick), HMC Step is the molecular dynam- 
ics timestep used (during replays, the step was halved) and lr ra j is the trajectory 
length used. 



of the 8 4 on our 12 core machines. 

4 Probability Distribution Function of Q 

First of all, we checked the behavior of the PDF of the operator Q outside the 
Aoki phase (point at p = 4.0 and K = 0.18) in the three volumes considered (4 4 , 
6 4 and 8 4 ). As figure |2] clearly indicates, the operator Q behaves as an intensive 
operator, as expected, and its PDF tends to a Dirac delta at the origin as the volume 
increases, showing that both parity and P' symmetries are realized in the vacuum 
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of the physical phase. 




Figure 2: PDF of operator Q outside the Aoki phase for three different volumes. 



Let's see what happens inside the Aoki phase (P = 2.0 and K = 0.25 and for 
this case we will use data from the 4 4 and 6 4 volumes (the 8 4 inside the Aoki phase 
was very expensive for us). We see in figure [3] and figure @] how the PDF behaves 
quite differently depending on the volume and on the sector. However, when we 
compute the final result taking into account the weight of each sector, we see in 
figure [5] how all these different plots converge to a single peak of constant width, 
in remarkable contrast with the results in the physical phase reported in figure 12 
Therefore we expect, according to our discussion in Section 2.4, that not only parity 
but also P' will be spontaneously broken inside the Aoki phase. We should notice 
that this result needs to hold if the standard picture of the Aoki phase is correct, 
since otherwise (cq) = (c|) / in this phase. 

Of course we are far from the thermodynamic limit, and one might argue that, 
as stated in Table [TJ the contribution of sector 2 to the 6 4 volume is important 
enough to be taken into account. This sector was not simulated because it was 
extremely expensive from the numerical point of view, and since its weight is less 
than 5%, we don't expect any changes to be relevant to the final result. 

Another indication that there are new vacua inside the Aoki phase not con- 
sidered in the standard picture, comes from measuring the PDF of the operator Q 
on configurations dynamically generated with several values of an external twisted 
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Figure 3: PDF of operator Q inside the Aoki phase for the two sectors considered 
in the 4 4 case. The sector 1 clearly shows a double peaked behavior 



mass term, m f /\j/y5X3\|/, but keeping m t = in the definition of Q. Since the addi- 
tion of an external source selects a standard Aoki vacuum, we expect that all the 
contributions to the PDF coming from the other vacua will be removed. As we can 
see in figure [6j the PDF for this case is essentially different than the one shown in 
figure [5] 

First of all, the PDF seems to be independent of the value of the external field. 
We might appreciate a slight dependence on the volume, since it seems that the 
peak height decreases as V increases, but this effect is too small to be significant, 
and in any case it is a good approximation to say that the distribution is also in- 
dependent of the volume. Comparing this distribution with the former one (Gibbs 
state, without external source, figure|5]), we notice that the PDF of the operator with 
external source is much wider, and we then expect the spectrum to be different as 
well. The essential difference lies in the low modes of the Dirac operator: in the 
case without external source, there is a lower bound given by 1 /V for the modu- 
lus of any eigenvalue, whereas at m, ^ the eigenvalues can become arbitrarily 
small H. Since the PDF depends strongly on the spectrum of the Dirac operator, we 

3 A consequence of this property is the fact that the PDF of Q at m, ^ fits perfectly to a 
Lorentzian distribution of infinite tails. One could also be concerned that Q become singular if a 
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Figure 4: PDF of operator Q inside the Aoki phase for the two sectors considered, 
now in the 6 4 volume. The double peak of sector 1 becomes narrower, but we 
should expect wider contributions from higher sectors. 



expect the cases m t = and m t / to be inherently different. 

We think that this result is the strongest indication of additional vacua structure 
in the Aoki phase outside the standard picture. Indeed the results of figure |7]clearly 
show that the sample of gauge configurations obtained in the Aoki phase at m t = 
is qualitatively different from the sample obtained at m t / 0, even in the m t — > 
limit. On the other hand the differences in the samples can never come from 
the additional standard Aoki vacua since if we change the twisted mass term in 
the dynamical generation by any other symmetry breaking term selecting other 
standard Aoki vacuum, for instance im^^^Xi^, the fermion determinant does not 
change. 



5 Second moments of the PDFs of fermion bilinears 

Since P' is broken according to our results plus the assumption that Q is an in- 
tensive operator, the symmetry arguments supporting the existence of a tower of 

zero-mode appeared in the spectrum. Actually, the set of configurations with a zero-mode has mea- 
sure zero, and Q remains non-singular. 
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-0.6 -0.4 



1/V Tr (y 5 D) 



Figure 5: PDF of operator Q inside the Aoki phase for two different volumes. The 
differences in the area can be explained by realizing that sector 2 is missing in the 
6 4 volume, which should account for a 5% of the area. 



sum rules are no longer valid, thus it is natural to explore if the expectation value 
^(/\Jry5\|/) 2 ^ will be non-zero in the Gibbs state inside the Aoki phase. This observ- 
able is extremely difficult to measure, nonetheless we managed to obtain sensible 
results from our data, which we show in Table [3] 

The fourth column in this table refers to ^(/\Jry5T3\|/) 2 ^, the landmark of the 
Aoki phase. As we see, it is clearly non-zero for all the volumes, confirming that 
our simulations lie within the Aoki phase. The second column shows the values for 
^('¥«Y5¥«) 2 y which is an order parameter for parity, but it takes into account just 
one flavor (which we labeled u). This quantity should be non-zero inside the parity- 
breaking Aoki phase, regardless of the discussion of the new vacua. Finally, the 

most important observable, the flavor singlet pseudoscalar ^(/vfrysXj/) 2 ^, which - 
since we are dealing with two degenerate flavors, u and d- is the sum of two of the 
former condensates rtjVY5Y« + ^WYsYd- The standard picture of the Aoki phase 
predicts zero expectation value of this parameter in any Aoki vacuum, whereas 
each one of the pseudoscalars restricted to one flavor Wfu4l5^u,d will be non zero 
due to parity breaking. Hence, the standard picture of the Aoki phase enforces 
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^4 
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A 
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• 


0.10 


A 
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▲ 
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0.10 


o 


0.30 


♦ 


0.50 





-0.5 

1/V Tr (Y 5 D)" 



Figure 6: PDF of operator Q inside the Aoki phase for configurations generated 
with a twisted mass source. 



an antiferromagnetic ordering of the pseudoscalars i\j/ M Y5¥« = ~~ z'WY5Vrf> which is 
not required in our hypothesis of the new vacua. The data of the third column in 
Table [3] show a clear non-zero expectation value in the case of the largest volume 
6 4 , supporting our previous discussion regarding P' breaking. 

Concerning equation (l22l ). which relates the mean value of the square Aoki 
condensate in the Gibbs state (s-regime) with the expectation value of this order 
parameter in the vacuum selected by a twisted-mass term (^-regime), and since we 
have data for the Aoki condensate at several lattice sizes and twisted masses (see 
Table @]), we can check its plausibility. This is a relevant test since, as discussed in 
Section 2.3, equation (1221 should be realized in the standard Aoki scenario 1221 12711 . 
To this end we report in figure [8]our results for the Aoki condensate at several val- 
ues of the twisted-mass, m t , in 4 4 ,6 4 and 8 4 lattices. The circles in the ordinates 
axis stand for the values of the Aoki condensate at m t = in 4 4 and 6 4 lattices, 
obtained from equation (l22l . and using as input our results for ((/\j/y5X3\j/) 2 ) re- 
ported in Table [3] The figure, which is a Fisher plot 14011 . shows that any reliable 
extrapolation of the data to m t = gives a value for the Aoki condensate larger 
than the one obtained from equation (l22l . 

One could argue that, at the volumes we are working, the finite size effects 
should be large and noticeable, and that the inclusion of these effects could lead to 
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Volume 


((i? a Y5¥«) 2 ) 


(('WsV) 2 ) 


((itj/ysTwf) 


<e 2 } 


4 4 Sec 
6 4 Sec 
4 4 Sec 1 
6 4 Sec 1 


(1.93 ±0.02) 10~ 2 
(2.15 ±0.05) 10~ 2 
(6.0±0.5)10- 3 
(1.75±0.09) 10~ 2 


(2.51 ±0.07) 10~ 2 
(2.82±0.18) 10~ 2 
(-5.5 ±0.5) 10~ 2 
(0.4 ±0.4) 10~ 2 


(5.23 ±0.09) 10~ 2 
(5.76±0.13) 10~ 2 
(7.9 ±0.4) 10~ 2 
(6.56±0.26) 10~ 2 


(6.8 ±0.3) lO" 3 
(7.4 ±0.6) ur 3 
(3.34±0.20) 10~ 2 
(1.53±0.15) 10~ 2 


4 4 Total 
6 4 Total 


(1.50 ±0.03) 10~ 2 
(1.81 ±0.12) 10~ 2 


(-0.6 ± 1.9) ur j 
(1.27±0.29) 10~ 2 


(6.08±0.19) 10~ 2 
(6.0±0.3) 10~ 2 


(1.54 ±0.08) 10" 3 
(1.17±0.13) 10~ 3 



Table 3: Interesting v.e.v. for the Aoki phase in the different sectors. Only statisti- 
cal eiTors are shown. 



V 


m t 


(l'W3Y) 


4 4 


0.01 


(56.00 ±0.06) lO -1 




4 4 


0.05 


(65. 149 ±0.022) 


10 


-2 


6 4 


0.01 


(53.603 ±0.015) 


10 


-2 


6 4 


0.05 


(66.888 ±0.008) 


10 


-2 


8 4 


0.05 


(66.953 ±0.007) 


10 


-2 


4 4 


0.10 


(75.5 12 ±0.009) 


10 


-2 


6 4 


0.10 


(76.291 ±0.004) 


10 


-2 


8 4 


0.10 


(76.307 ±0.003) 


10 


-2 


4 4 


0.30 


(92.214 ±0.005) 


10 


-2 


6 4 


0.30 


(92.309 ±0.002) 


10 


-2 


8 4 


0.30 


(92.311 ±0.002) 


10 


-2 


4 4 


0.50 


(96.225 ±0.003) 


10 


-2 


6 4 


0.50 


(96.241 ±0.001) 


10 


-2 


8 4 


0.50 


(96.243 ±0.001) 


10 


-2 



Table 4: Evolution of the v.e.v. of the Aoki parameter as a function of the twisted 
mass external field. 
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Figure 7: PDF of operator Q inside the Aoki phase for all configurations. Compar- 
ison of the case m t ^ versus the case m t = 0. 



an agreement between our measurements at m t = and m t ^ 0. Nonetheless, the 
data at m t = -which we expect to suffer more these finite volume effects due to 
the existence of massless pions- reveals that these effects are not very large, for 
the value of the Aoki condensate remains very stable after a fivefold increase in the 
volume. 



6 Conclusions and Outlook 

Three years ago the standard scenario for the phase structure of lattice QCD with 
two degenerate flavors of Wilson fermions was revised by three of us in 12111 . where 
we conjectured on the appearance of new vacua in the Aoki phase, which can be 
characterized by a non-vanishing vacuum expectation value of ifyulstyu + i^ifl5^fd, 
and which can not be connected with the Aoki vacua by parity-flavor symmetry 



transformations. However, Sharpe pointed out in [22] that the standard picture for 
the Aoki phase can be understood using the chiral effective theory appropriate to 
the Symanzik effective action, and that within this standard analysis, the flavor- 
singlet pseudoscalar expectation value vanishes, (i\j/y5\|/) = 0. As the standard 
scenario for the Aoki phase is indeed a direct consequence of the %PT application 
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Figure 8: Fisher plot of the evolution of the v.e.v. of the Aoki parameter as a 
function of the twisted mass external field. The circles show the prediction obtained 
from the simulations in the Gibbs state (see eq. (l22l) ). The fit shown is meant to 
guide the eye, and shouldn't be taken as a prediction of the behavior of the v.e.v. 
of the Aoki parameter. The fit should be more reliable for small values of m t . 



to this phase, we were also calling into question in j2lll the validity of the %PT 
analysis, at least for low values of p\ and therefore large values of the lattice spac- 
ing a. These issues are relevant enough to make it worthwhile to continue these 
investigations in order to clarify the actual scenario for the Aoki phase. 

For the last few years we have performed an extensive research on the vacuum 
structure of the Aoki phase, in order to find out if our alternative vacuum structure, 
derived from a PDF analysis, was realized or not. These simulations, which have 
been mainly performed in the absence of any parity-flavor symmetry breaking ex- 
ternal source, are plagued by technical difficulties which have been responsible for 
the slow progress in the field. Indeed the addition of a twisted mass external source, 
as it has been done in the past simulations of the Aoki phase, automatically selects 
the vacuum where the standard properties of the Aoki phase are verified. This 
point could explain why nobody ever saw a new phase like the one we are propos- 
ing, since all the past simulations performed with two flavors of Wilson fermions 
inside the Aoki phase were done within an external twisted mass term. 

Notwithstanding these difficulties, we have provided in this work evidence 
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pointing to a more complex vacuum structure in the Aoki phase of two-flavor QCD, 
as conjectured in |2lll . Indeed the results reported in Table [3j which were obtained 
under the assumption that the weights of the different sectors are continuous at 
m t = 0, show how we were able to perform a direct measurement of ^(/\j/y5\|/) 2 ^) 
in the Aoki phase, which gave a non-zero value for this operator in the 6 4 lattice, 

of the same order of magnitude as (^(iyufstyu)/, an d ^('W5 X 3V) 2 ^)> tne l ast two 
being non-vanishing in the Aoki phase in the standard scenario. Thus this result 
points to the breaking of the hypothesis of the sum-rules (fT3l) . Furthermore a check 
of equation (1221) . a equation which should hold if the standard scenario for the Aoki 
phase is realized, reported in figure 8, points out to a more complex vacuum struc- 
ture too. 

However the strongest indication, in our opinion, on the existence of a vacuum 
structure in the Aoki phase, more complex than the one of the standard picture, 
comes from our analysis of the PDF of the operator Q (1231) . Our motivation for the 
analysis of this kind of density of "topological charge" operator, which measures 
the asymmetries in the eigenvalue distribution of the hermitian Dirac -Wilson op- 
erator, was twofold. First it appears in the computation of the second moment of 
the PDF of the flavor singlet pseudoscalar order parameter (fT2l . and second Q is 
an order parameter for the symmetry P', composition of parity and discrete flavor 
transformations, described in Section 2. 

Our numerical results for the PDF of Q, together the assumption that Q is an 
intensive operator, suggest that the P' symmetry is realized in the vacuum of the 
physical phase, but not in the Aoki phase, at least in the strong coupling region 
we have analyzed. However the cleanest signal of further structure in the Aoki 
phase comes from the results on the PDF of Q reported in figure [7] Those results 
clearly show that the sample of gauge configurations obtained in the Aoki phase 
at m t = is qualitatively different from the sample obtained at m t ^ 0, since they 
give incompatible PDF's. Furthermore this result stays true in the m t — > limit, 
as follows from the independence of the shape of the PDF of Q on the twisted 
mass m t . On the other hand the differences in the samples can never come from 
the additional standard Aoki vacua since if we change the twisted mass term in 
the dynamical generation by any other symmetry breaking term selecting other 
standard Aoki vacuum, as for instance imt^ysXi^, the fermion determinant does 
not change. 

As the sum rules (fT3l) required in the standard scenario are only supported, as 
discussed before, by the %PT analysis of the Aoki phase, our results could call into 
question also the validity of the %PT analysis performed in HI lfl for low values of 
P (around 2.0, very coarse lattices). However, %PT is expected to work at higher 
values of p\ near the continuum limit. Indeed, one could be concerned that, at such 
a low p\ we are far from the continuum limit. Nonetheless our analysis is devoted 
to the Aoki phase on the lattice, which might not even have a continuum limit. 

In any case any improvement of our results in larger lattices and for more (3,K 
values would be of great interest. Unfortunately, given the technical difficulties 
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of the simulations inside the Aoki phase, this calculation is outside our present 
computing resources. 
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